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DISCLAIMER 



Although each program has been tested by its contributor, no 
warranty, express or implied, is made by the contributor or 
1620 USERS Group, as to the accuracy and functioning of the 
program and related program material, nor shall the fact of ^ 
distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or 1620 USERS Group, in con- 
nection therewith. 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 
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Commen t 

Are the Sense Switch options adequately described (if applicable)? Yes No 

Are the mnemonic labels identified or sufficiently understandable? Yes No 

Commen t 

5. Does the source program compile satisfactorily (if applicable)? Yes No 

Commen t 

6. Does the object program run satisfactorily? Yes _No 

Comment 

7. Number of test cases run . Are any restrictions as to data, 
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Commen t 
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ABSTRACT 
LINEAR PROGRAMMING II 
by 

F. H. flood 



DESCRIPTION 



This Linear Programming code is essentially the same as the previous 
four-phase code, the difference being the addition of an output phase and 
the reduction in the number of typeouts in the other phases. This writeup 
assumes familiarity with the four-phase writeup, the difference being what 
is included here. 



DISCUSSION 



The revision of the four-phase program to the five-phase program was 
initiated by the need for extra matrix capacity. Error typeouts were 
eliminated and only the pause now exists. The output portion of the program, 
after the optimum solution is reached, was removed from the Solution Phase 
and from the Right Hand Side Changer Phase and made an Individual Phase. 
Therefore, to obtain a solution it is necessary to load the Output Phase after 
the final stop has been encountered in the Right Hand Side Chdnger or the 
Solution Phase. 

The MAR display number as described in the following sections are all 
based on FORTRAN with auto-divide and a lowered origin. If the programs are 
compiled with FORTRAN decks that do not comply with the above, these numbers 
will differ from those reported below and the correct number must be established 
from the print off when the program is compiled. 

DATA LOADER PHASE 

Input 

Same as L.P.I. 
Output 

Same as L.P.I. 
Sense Switch Settings 

On Off 
SS 1 For restart problems For initial start problems 

SS 2,3,4 Not used 
Error Stops 
None 



NATIONAL STEfcl RESEARCH 



r 



2 - 



SOLUTION PHASE 
Input 

None 
Output 

Same as L.P, 1 except no headings and no final output. 
Sense Switch Settings 

On Off 



SS 4 



SS 1,2,3 
Error Stops 
MAR Display 
9039 
12011 
12023 
COST CHANGER PHASE 
Input 

Same as L.P. 1 
Output 

Same as L.P. 1 
Sense Switch Settings 

None used 
Error Stops 

None 

RIGHT HAND SIDE CHANGER PHASE 
Input 

Same as L.P. 1 



Punch output for restart at a 
later time 

Not used 



Continue with Simplex 
Method 



Meaning 

Problem as stated is unbounded 
Problem as stated is infeasible 
Optimum solution has been reached 



NATIONAL 



STEEL 



R E S E. A R 
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RIGHT HAND SIDE CHANGER PHASE (Cont'd) 
Output 

Same as L.P. 1 with the exception of error typeouts and the 
final solution. 

Sense Switch Settings 

None used 
Error Stops 

MAR Display Meaning 

10651 Problem is infeasible 

12575 Optimum feasible solution has been reached 



An alternate method of input is available without any program change. 
The original method of input was to punch both the old right hand side element 
and the new right hand side element into the change card. 

It is possible to record only the difference betweenthe elements in card 
columns 9-18 and leave the rest of the card blank. 



The card should then be punched as follows: 



cc 

1-4 l) 0000 if the unit vector associated with tne change is artificial 

2) 0001 if the unit vector associated with the change is slack 

5-8 Column number of unit vector associated with the change 

9-18 Old right hand minus the new right hand side with a decimal 

point somewhere in the 10 digit field. Care must be taken 
with the sign of the resultant of the subtraction. To 
increase the right hand side element, the number must be 
negative and to decrease it must be positive. 

OUTPUT PHASE 
Input 



Output 

Solutions shadow costs, final matrix punch, and solution punch 



NATIONAL STEEL RESEARCH 
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f 



r 



Sense Switch Settings 
On 

SS 1 Print final solution 

SS 2 Punch final matrix 

SS 3 Print Shadow Costs 

SS 4 Punch final solution 

In all the above cases the final 
Error Stops 
None 



Off 

Do not print final solution 
Do not punch final matrix 
Do not print Shadow Costs 
Do not punch final solution 

solution is printed on the typewriter. 



Mathematician 
FWW/jd 



O N A L 



STEEL 
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DAT. 



-OADER 



:(;CO 
C . 00 
G.'OCC 

0-023 

OwOCO 
C'M 14 
£■150 
C3i36 
c:i : j2 
00246 

03294 
or; 3 06 
0H462 

g n 5 1 o 

03^66 
08350 
Jb 5 o 
03632 
03730 
03736 
C Li 8 70 
03913 
03926 
O0-;4S 
03970 
09006 
C9074 
0515^ 
091 66 
09282 
09493 
C953+ 
09542 
09566 
09578 



D i MEi'S I ON A(44,99) , ' / ( 44 ) , L ( 44 ) , I I ( 1 ) .JJ ( 1 ) , I I I ( 1 ) ,X( 1 ) 
DIMENSION JK(1),KJ(1),K:<K(1),IJ(1),k(i),I(1),J(1),XMIN(1) 

1 F0RMAT(l4,l4) 

2 FORMAT ( I 4 , 1 4 ,F20 .8) 

4 FORMAT (15HL0AD NEXT FHASE) 
108 PEAD 1,11 , JJ 
111=11+1 
DO 10 1=1 ,1 I I 

V/( I )=0.0 - - 

L(l)=0 

DO 10 J = 1 ,JJ 

10 A( I ,J)=0.0 

11 READ 2, 1 ,J,X - 
I F ( I ) 16,16,18 

18 A( I ,J)=X . - 

GO TO 11 

16 J=JJ - - -■ - 

17 READ 2, I ,JK,X 

I F ( I ) 20 , 20 , 1 9 - - — 

19 A( I ,JJ)=X 

L(I)=JK — - - 

GO TO 17 

20 I F ( SENSE SWITCH 1)40,22 

22 1=1 

23 1=1+1 " " " - - 
IF( l-l I I ) 24 ,40 ,40 

24 IF(L(! ))23,25 ,23 

25 DO 27 J = 1 ,JJ 

IF(A(I ,J))26,27,26 

26 A( I I I ,J)=A( I I I ,J)-A( I ,J) 

27 CONTI NUE -■ ■ • 

GO TO 23 

40 PRINT 4 -- —- 

PAUSE 

END - - 
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SOLUTION PH 



03000 
03 COO 
03000 
03032 
03056 
3030 
03104 
03152 
03200 
03236 
03304 
08420 
03572 
03600 
03723 
08736 
03816 
08864 
08876 
03992 
09028 
09040 
09064 
09038 
09124 
09192 
09308 
09464 
09520 
09588 
09612 
09636 
09644 
09728 
09776 
09788 
Q9932 
09968 
09980 
09992 
10108 
10188 
10464 
10536 
10572 
:-o34 
10596 
10712 
10792 
11068 



I MENS I OH A(44,99) ,1/(44) , 1.(44), I I (1),JJ(1),I I I ( 1 ) , X ( 1 ) 

D I ME MSI ON JK( 1 ) ,KJ ( 1 ) ,KKK( 1 ) , IJ (1 ) ,K(1 ) , I ( 1 ) , J ( 1 ) ,XMI N( 1 ) 
105 FORMAT (I4,F15.2,I4) 

KKK=0 

i\= I I I - 

44 J=0 
W(K)=0.0 
L(K)=0 

42 J=J+1 " 
!F(J-JJ)41 ,45,45 

41 IF(A(!<,J))43,42,42 

43 IF(W(K)-A(K,J))42,42,47 
47 >.'(K)=A(K,J) 

L(K)=J 
GO TO 42 

45 IF(L(K))46, 62,46 

46 KJ=L(K) 

DO 120 1=2,1 I 

IF(A( I ,KJ))120,120,121 - - ~ - 

120 CONTINUE 
PAUSE 

121 1=1 

JK=0 - 

50 1=1+1 

1 F( I - 1 I )52,52,56 

52 I F( A( i ,KJ)-)50, 50,51 

51 X=A( I ,JJ)/A(I ,KJ) 

IF(JK)55,53,55 

55 IF(X-XMIN)53,50,50 

53 XMIM=X 

GO TO 50 

56 X=A(JK,KJ) - - 

L(JK)=KJ 

DO 57 1=1,111 

57 V/(l )=A(I ,KJ) 
I J=JK-1 

DO 59 1=1 , IJ 
DO 59 J = 1 ,JJ 
IF(A(JK,J))58,59,58 

58 I F (V/( I ) ) 580 ,59 ,580 

580 A(l ,J)=A(I ,J)-W(I)*(A(JK,J)/X) 

59 CONTINUE 
IJ=JK+1 

DO 61 1 = 1 J, I I I 
DO 61 J=1 ,JJ 
IF(A(JK,J) )60,61 ,60 

60 IF(W( I ))600,61 ,600 

600 A(l ,J)=A(I ,J)-W(I)*(A(JK,J)/X) 

61 CONTINUE 



SOLUTION PHAS-: (CGi-!T) 



JM 140 DO 205 J=1 ,JJ 

11152 205 A(JK,J)=A(JK,J)/X 

7 1344 KKK=KKK+1 

11380 PRINT 105,KKK,A(K,JJ) ,L(JK) 

11512 IF(SENSE SWITCH 4)70,44 

11532 62 IF(K-1)70,70,63 

11 600 63 IJ=JJ-1 

11636 DO 65 J=1 ,IJ 

7l648 IF(A(K,J)-.0001)65,65,66 

11776 65 CONTINUE 

11312 DO 130 J=1,JJ 

T1B24 130 A( I I I ,J)=0.0 

11944 K=1 

11963 KKK=0 

11392 GO TO 44 

12000 66 PAUSE 

12012 70 PAUSE 

12024 END 



e 



IGHT HAHD SIDE CHANGER 



c;:ooo dikemsiom a(44,99) ,w(44) ,l<44) , i 1 d ) ) , 1 1 1 (1 ) ,x(i ) 

OoOOO D I MEi IS I ON JK{ 1 ) ,KJ ( 1 ) ,iC;K( 1 ) , U ( 1 ) ,K( 1 ) , I ( 1 ) , J ( 1 ) ,XMI N( 1 ) 

COGOO 1 FORMAT (I4,I4,F10.5,F10.5)' 

08033 3 FORMAT (33H ITER FUNCTIONAL VAR. IN) 

■03123 100 F0RMAT(F20.5,F20.5) 

00156 101 FORMAT (32H OLD PHS NEW RHS) 

08244 1 04 FORMAT ( i 4 , F 1 5 • 2 , 1 0X , I 4 ) 

03304 8 PRINT 101 

03328 K=1 ' ' 

00352 KKK=0 

08376 4 READ 1 ,KJ , JK ,W( 1 ) ,W( III) ~ ~ 

08460 IF(KJ)01 ,80,81 

08516 80 IF(JK)82,11 ,82 - 

08572 81 IJ=1 

: 08596 GO TO 5 " 

03604 82 IJ=2 

08628 5 PRINT 100 ,V/( 1 ) ,W( III; 

08688 V/(1)=W(M l-)-W(l) 

03748 DO 10 l=IJ,l I " • 

68760 10 A(l ,JJ)=A(I ,JJ)+W(1)*A(I,JK) 

09024 GO TO 4 ' 

09032 11 PRINT 3 

09056 15 XMIN=0.0 '-- ~ 

09080 L(K)=0 

09128 1=1 ~ 

09152 18 [=1+1 

09188" IF( f-1 1)16,16,20 — — "-• — — • — 

09256 16 IF(A(t,JJ))17,fS,l8 

09372 17~!F(A(1 ,JJ)-XMIN)19, 18,18 ~— ~~ . - . 

09500 19 XMIM-A(I ,JJ) 

09584 - L(K) = I • 

09632 GO TO 18 

09640 20 IF(L(K))21 ,70,21 " ' 

09720 21 JK=L(K) 

09768 ~ " J=0 

09792 XMIM=0.0 

09816 " L(K)=0 " 

59364 26 J=J+1 

09900 IF(J-JJ)23 ,24,24 " " ' — 

09968 23 IF(A(JK,J))25,26,26 

70084 25 IF(A(1 ,J))26,26,7 ' ' ' 

T0188 7 X=A(1 ,J)/A(JK,J) 

10332 IF(L(K))6,27,6 ' • 

70412 6 1F(X-XMIH)26,26,27 

10480 '" " 27 XK.IN=X ' " 

T0504 L(K)=J 

"10552 " ' GO TO 26 

10560 24 IF(L(K))31 ,30,31 

10640 30 PAUSE ' ~: • 

70652 31 KJ=L(K) 



JO 



o 
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RIGHT HAND Si DE CHANGER (CONT) 



10700 L(JK)=KJ 

70748 X=A(JK,KJ) 

1C832 DO 33 1 = 1 ,1 I I 

TO 344 33 '/(I )=A( I ,KJ) 

T0908 IJ»JK-1 

11024 DO 35 1=1 ,U 

J[1036 DO 35 J=1 ,JJ 

11048 IF(A(JK,J))34,35,34 

31164 34 IF(W(I))37,35,37 

11244 37 A(l ,J)=A(I ,J)-W(I)*(A(JK,J)/X) 

11520 35 CONTINUE 

11592 IJ=JK+1 

T1628 DO 40 1 = 1 J, I II 

11S40 DO 40 J-1.JJ 

T1652 IF(A(JK,J))38,40,38 

T1763 38 IF(W(I ))39,40,39 

11843 39 A(l ,J)=A(I ,J)-W(l)*{A(JK,J)/xr 

T21 24 40 CONTINUE 

12196 DO 50 J=1 ,JJ 

12208 50 A(JK,J)=A(JK,J)/X 
12400 KKK=KKK+1 

T2436 PRINT 104,KKK,A(1,JJ) ,L(JK) 

T2556 GO TO 15 

12564 70 PAUSE 

T2576 GO TO 8 " 

12584 END 
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COST CHANCER 



onooo 
gaooo 

03000 
03032 
0G082 
OB132 
081 92 
08266 
08290 
08362 
O8386 
08442 
08514 
O8574 
08706 
08718 
08810 
68846 
08854 
08866 
08982 
09198 
09234 
09242 
09278 
09290 
Q9394 
09430 
69454 
09526 
09538 
09546 
09570 
09582 
09590 



D I MEMS I ON A(44,99) .W( 44) ,L(44) , I I ( 1 ) ,JJ( 1 ) , I I I ( 1 ) ,X( 1 ) 
D I MENS I ON JK( 1 ) ,KJ ( 1 : ,KKK( 1 ) , IJ ( 1 ) ,K( 1 ) , I ( 1 ) , J ( 1 ) ,XMI N( 1 ) 

1 FORMAT ( I4,F10.5»F10.5) 

2 FORMAT (13HST ILL OPTIMAL) 

3 F0RMAT(10HFUNCTI0NAL,F20.8) 

4 FORMAT (18HL0AD SOLUTION DECK) 

100 FORMAT (25HVAR. OLD COST NEW COST ) 
9 PRINT 100 

10 READ 1 ,JK,W(1),W( III) " 
L(1)=JK 

I F( JK) 11 ,50 ,1 1 

11 PRINT 1 ,JK,W(1),W( III) 
V/(1)=W( I I I ) -W ( 1 ) 

A(1 ,JK)=A(1 ,JK)+W(1) 

DO 15 1=2,1 I ' • ""' 

IF(L(I)-L(1))15.16,15 

15 CONTINUE : - " "' 
GO TO 10 

16 DO 20 J=1 ,JJ " ""' " 
IF(A(I ,J))19,20,19 

19 A(1,J)=A(1,J)-W(1)*A(I,J) ' ' " 

20 CONTINUE 

GO TO 10 " 
50 I J=JJ-1 
, - DO 55 J=1 ,U 

IF(A(1,J))60,55,55 

CONTINUE ~ 



55 



60 



PRINT 
PRINT 
PAUSE 
GO TO 
PRINT 
PAUSE 
GO TO 
END" 



3,A(1 ,JJ) 

9 
4 



r 



OUTPUT ' J 



g.iooo 

g:<CCO 

g3ogo 

C1C32 

g3060 
g:;i30 

G3250 

03276 

3296 
C .I360 
03332 

0"o60 
03530 
03604 
OGohO 
O0&52 
OC-756 
0.3840 
03376 
g833& 
03932 
08944 
63956 
03072 
031 80 
09252 
03276 
09288 
09^56 
03468 
09488 



e 



D I ME 3 S I ; ! A ( 44 , 9 9 ) , ' •/ \ 44 ) , L ( 44 ) , I I ( 1 ) . J J ( 1 ) , I i I ( 1 ) , X ( 1 ) 
D I MEMS I OH JK(1) ,KJ ( 1 ) ,KKK( 1 ) , IJ (1 ) ,K( 1) , I ( 1 ) , J ( 1 ) ,XMI N( 1 ) 

2 FORMAT ( 14, I4.F20.8) 

5 FORMAT ( I4.F20.3) 

7 FORMAT ( 23HVARI ABLE VALUE) 

8 FORMAT ( 1 OHFUNCT I ONAL , F20.3) 

9 FORMAT ( 23 HVAR I ABLE SHAD. COST) 
100 FORMAT ( 1 H ) 

I F( SENSE SWITCH 1)70,92 

70 PR I NT 8,A( 1 ,JJ) 
PRINT 7 

DO 71 1=2, 1 I 

71 PRINT 5,L(I),A(I ,JJ) 

92 I F( SENSE SWITCH 3)72,75 

72 PRINT 9 
! J=JJ-1 

DO 74 J = 1 ,IJ 
IF(A(1 ,J))73,74.73 

73 PRINT 5, J ,A(1 ,J) 

74 CONTINUE 

75 I F ( SENSE SWITCH 2)76,80 

76 I J=J J-1 

DO 77 J = 1 ,U 
DO 77 1=1 J i i 
IF(A(I ,J))78, 77,78 
78 PUNCH 2,1 ,J,A(I ,J) 

77 CONTINUE 
PUNCH 100 

91 DO 90 1 = 1,111 " ~ 

90 PUNCH 2, I ,L( I ) ,A( i ,JJ) 
81 PAUSE 

80 IF(SENSE SWITCH 4)91 ,81 
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